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This article describes a determination of the Cabibbo-Kobayashi-Maskawa matrix element \V c b\ 
from the decay B° — > D*~ l + vi using 711 fb _1 of Belle data collected near the T(45*) resonance. We 
simultaneously measure the product of the form factor normalization J-(l) and the matrix element 
1 14& | as well as the three parameters p 2 , i?i(l) and i?2(l), which determine the form factors of this 
decay in the framework of the Heavy Quark Effective Theory. The results, based on about 120,000 
reconstructed B° ->■ D*~l + vi decays, are p 2 = 1.214 ± 0.034 ± 0.009, i?i(l) = 1.401 ± 0.034 ± 0.018, 
Jfe(l) = 0.864 ± 0.024 ± 0.008 and J^(l) | V c& | = (34.6 ± 0.2 ± 1.0) x 10" 3 . The branching fraction 
of B° — > D*~l + ve is measured at the same time; we obtain a value of B(B° —>■ D*~ l + Vf) = 
(4.58 ± 0.03 ± 0.26)%. The errors correspond to the statistical and systematic uncertainties. These 
results give the most precise determination of the form factor parameters and .F(l)|V c b| to date. 
In addition, a direct, model-independent determination of the form factor shapes has been carried 
out. 



I. INTRODUCTION 

The study of the decay B° — > D*~£ + ve is important 
for several reasons. The total rate is proportional to the 
magnitude of the Cabibbo-Kobayashi-Maskawa(CKM) 
matrix element V c b [H, squared. Experimental investi- 
gation of the form factors of the decay can check theoret- 
ical models and possibly provide input to more detailed 
theoretical approaches. In addition, B° — > D*~£ + ve is a 
major background for charmless semileptonic B decays, 
such as B — > -nlv, or semileptonic B decays with large 
missing energy, including B — > D*tv. Precise knowledge 
of the form factors in the B° —> D*~l + vi decay will thus 
help to reduce systematic uncertainties in these analyses. 

This article is organized as follows: After introduc- 
ing the theoretical framework for the study of B° — > 
D*~£ + vg decays in Section UH the experimental proce- 
dure is presented in detail in Section IIIII This is fol- 
lowed by a discussion of our results and the systematic 
uncertainties assuming the form factor parameterization 
of Caprini et al. [1| in Section IIV1 Finally, a measure- 
ment of the form factor shapes is described in Section fVl 

This paper supersedes our previous result (36j . based 
on a subset of the data used in this analysis. 



II. THEORETICAL FRAMEWORK 
A. Kinematic variables 

The decay B° -t D*~l + vt Q proceeds through the 
tree-level transition shown in Fig.Q] Below we will follow 
the formulation proposed in reviews [f| @, where the 
kinematics of this process are fully characterized by four 
variables as discussed below. 

The first is a function of the momenta of the B and 




FIG. 1: Quark-level Feynman diagram for the decay B 
D*~i + u e . 



D* mesons, labeled w and defined by 



w 



D * 



niBrriD* 



2m,BmD' 



(1) 



where tub and mr)* are the masses of the B and the D* 
mesons (5.279 and 2.010 GeV/c 2 , respectively [?]), P B 
and Pd* are their four-momenta, and q 2 = (Pi + P v ) 2 . 
In the B rest frame the expression for w reduces to the 
Lorentz boost 7d» = Ed* jmu* ■ The ranges of w and 
q 2 are restricted by the kinematics of the decay, with 
q 2 ■ 

^min 



corresponding to 

™2 



277JbTO£|. 



and 



= 1 to 



(m B - m D * Y 



1.504 , 



10.69 GeV 2 



(2) 



(3) 



The point w = 1 is also referred to as zero recoil. 

The remaining three variables are the angles shown in 

Fig.m 

• 6 1, the angle between the direction of the lepton 
and the direction opposite the B meson in the vir- 
tual W rest frame; 
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FIG. 2: Definition of the angles 61, 9v and \ f° r the decay 
B° — > D*~ t + v f , D*~ — D°tv7 . 



• 6>y , the angle between the direction of the D meson 
and the direction opposite the B meson in the D* 
rest frame; 

• Xi the angle between the plane formed by the D* 
decay and the plane formed by the W decay, de- 
fined in the B meson rest frame. 



B. Four-dimensional decay distribution 

Three helicity amplitudes, labeled H + , H-, and Ho, 
can be used to describe the Lorentz structure of the 
B° —> D*~l + vi decay amplitude. These quantities cor- 
respond to the three polarization states of the D* , two 
transverse and one longitudinal. When neglecting the 
lcpton mass, i.e., considering only electrons and muons, 
these amplitudes are expressed in terms of the three func- 
tions hA 1 (w), Ri(w), and R 2 (w) as follows @: 

R*(l-r 2 )(w + l) u - 

Hi(w) = m B - r h Al (w)Hi(w) , (4) 

2y 1 — 2wr + H 

where 



Vl-2wr + r* (i± y /%g Rl ( w ] 



Ho = 1 



1 -T 

(w - - R 2 (w)) 
1 - r 



(5) 
(6) 



with R* = (2^/msm£)*)/( m s+TO_D») and r = itid* /tub- 
The functions R\(w) and R 2 (w) are defined in terms of 
the axial and vector form factors as, 



A 2 (w) 



V{w) 



R 2 (w) 



>*2 



R 



Ri(w) 
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M{w) . 



(7) 



(8) 



By convention, the function h Al (w) is denned as 

1 2 



h Al (w) 



R* w + 1 



Ax(w) 



(9) 



The axial form factor A\(w) dominates for w 1. Fur- 
thermore, in the limit of infinite &- and c-quark masses, 
a single form factor describes the decay, the so-called 
Isgur- Wise function d, HI ■ 

In terms of the three helicity amplitudes, the fully dif- 
ferential decay rate is given by 



d 4 T(B° -> D* 



~vi) 



6mBm 



dw d(cos6 e ) d(cos9 v ) dx 8(4tt) 
(1 - 2wr + r 2 )G 2 F \V cb \ 2 {(l -cos6» £ ) 2 sin 2 6 v H%{w) 
+ (1 + cos 9 e ) 2 sin 2 9 V H 2 _ (w) + 4 sin 2 9 e cos 2 9 v H 2 (w) 

- 2 sin 2 ^ sin 2 6V cos 2xH + (w)H- (w) 

— 4 sin 6*^(1 — cos^) sin#y cos6>y cos xH+(w)Ho(w) 

+ 4sin $t(l + cos9e) s'm9v cos9v cos xH-(w) Ho (w)} , 

(10) 

with G F = (1.16637 ±0.00001) x 10~ 5 ?ic 2 GeV~ 2 . Four 
one-dimensional decay distributions can be obtained by 
integrating this decay rate over all but one of the four 
variables, w, cos 9g, cos9y, or x- The differential decay 
rate as a function of w is 



dr _ Gj 

dw 487i 
where 

F 2 {w)g{w 



mLK - m D ,) 2 g{w)F 2 {w)\V cb \ 2 , (11) 
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l + i?i(w) 5 



w - 



1 + (1-R 2 (w)) 



l-r 



and Q{w) is a known phase space factor, 

w 1 — 2wr + r 2 1 



G(w) = \Jw 2 - l(w + if 



1 + 4 



w + 1 (1 - r) 2 



A value of the form factor normalization .FQJ = 1 is 
predicted by Heavy Quark Symmetry (HQS) |5| in the 
infinite quark-mass limit. Lattice QCD can be utilized 
to calculate corrections to this limit. The most recent 
result obtained in unquenched lattice QCD is .F(l) 
0.921 ±0.013 ±0.020 



i que 
.10] 



C. Form factor parameterization 

A parameterization of form factors h Al {w) , R\ (w) , and 
R 2 (w) can be obtained using heavy quark effective the- 
ory (HQET). Perfect heavy quark symmetry implies that 
Ri(w) = R 2 (w) = 1, i.e., the form factors A 2 and V are 
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identical for all values of w and differ from A\ only by a 
simple kinematic factor. Corrections to this approxima- 
tion have been calculated in powers of Aqcd / "it and the 
strong coupling constant a s . Various parameterizations 
in powers of (w — 1) have been proposed. We adopt the 
following expressions derived by Caprini, Lellouch and 
Neubert 0], 

h Al (w) = h Al {l)[l~8p 2 z + (53p 2 -15)z 2 

- (231p 2 -91)z 3 ] , (12) 
Ri(w) = i?i(l)-0.12(w-l) + 0.05(w;-l) 2 ,(13) 
Ri{w) = i? 2 (l)+0.11(w-l)-0.06(w;-l) 2 ,(14) 

where z — (y/w + l — y/2)/(\/w + l + \/2). In addition to 
the form factor normalization .F(l) = hj^ x {V), these ex- 
pressions contain three free parameters, p 2 , i?i(l), and 
i? 2 (l). The values of these parameters cannot be calcu- 
lated in a model-independent manner. Instead, they have 
to be extracted by an analysis of experimental data. 



four distinct sets labeled A (141 ftr 1 ), B (274 ftr 1 ), C 
(189 ftr 1 ) and D (107 ftr 1 ), where the number in paren- 
theses indicates the integrated luminosity corresponding 
to the individual samples. 

Monte Carlo generated samples of BB decays equiva- 
lent to about three times the integrated luminosity are 
used in this analysis. Monte Carlo simulated events are 
generated with the Evtgen program [Til ], and full detec- 
tor simulation based on GEANT [15| is applied. QED final 
state radiation m_B -> Xlv decays is added using the 
PHOTOS package [l6f . 

Hadronic events are selected based on the charged 
track multiplicity and the visible energy in the calorime- 
ter. The selection is described in detail elsewhere [l7j . 
We also apply a requirement on the ratio of the second 
to the zeroth Fox-Wolfram moment 18], i?2 < 0.4, to 
reject continuum events. 



B. Event reconstruction 



III. EXPERIMENTAL PROCEDURE 
A. Data sample and event selection 

The data used in this analysis were taken with the 
Belle detector [ill at the KEKB asymmetric-energy 
e + e~ collider [liq . The Belle detector is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50- layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel-like arrangement of time-of- flight scintil- 
lation counters (TOF) , and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect mesons and to iden- 
tify unions (KLM) . The detector is described in detail in 
Ref. [ll|. Two inner detector configurations were used. 
A 2.0 cm beampipe and a 3-layer silicon vertex detector 
were used for the first sample of 152 million BB pairs, 
while a 1.5 cm beampipe, a 4-layer silicon detector and 
a small-cell inner drift chamber were used to record the 
remaining 620 million BB pairs [l3j |. 

The data sample consists of 711 ftr 1 taken at the 
T(45) resonance, or about 772 million BB events. An- 
other 88 ftr 1 taken at 60 MeV below the resonance are 
used to estimate the non-BB (continuum) background. 
The off-resonance data is scaled by the integrated on- to 
off-resonance luminosity ratio corrected for the 1/s de- 
pendence of the e + e~ — ¥ qq cross section. 

This data sample contains events recorded with two 
different detector setups as well as two different tracking 
algorithms and large differences in the input files used 
for Monte Carlo generation. To ensure that no system- 
atic uncertainty appears due to inadequate consideration 
of these differences, we separate the data sample into 



Charged tracks are required to originate from the in- 
teraction point by applying the following selections on 
the impact parameters in the r — <fi and z directions: 
dr < 2 cm and \dz\ < 4 cm, respectively. In addition, we 
demand at least one associated hit in the SVD detector. 
For pion and kaon candidates, the Cherenkov light yield 
from the ACC, the time-of-flight information from TOF, 
and dE/dx from the CDC are required to be consistent 
with the appropriate mass hypothesis. 

Neutral D meson candidates are reconstructed in the 
D° — > K~ir + decay channel. We fit the charged tracks to 
a common vertex and reject the D° candidate if the x 2 - 
probability is below 10~ 3 . The reconstructed D° mass is 
required to lie within ±13.75 MeV/c 2 of the nominal D° 
mass of 1.865 GeV/c 2 j7}, corresponding to about 2.5 
times the experimental resolution measured from data. 

The D° candidate is combined with an additional 
charged pion (oppositely charged with respect to the kaon 
candidate) to form a D* + candidate. Due to the kine- 
matics of the D* + decay, the momentum of this pion 
does not exceed 350 MeV/c. It is therefore referred to as 
the "slow" pion, 7r+. No impact parameter or SVD hit 
requirements are applied for ir s . Again, a vertex fit is per- 
formed and the same vertex requirement is applied. The 
invariant mass difference between the D* and the D can- 
didates, Am = mr>* —rnrjo , is required to be less than 165 
MeV/c 2 . This selection is tightened after the background 
estimation described below. Additional continuum sup- 
pression is achieved by requiring that the D* momentum 
in the cm. frame be below 2.45 GeV/c. 

Finally, the D* candidate is combined with an oppo- 
sitely charged lepton (electron or muon). Electron candi- 
dates are identified using the ratio of the energy detected 
in the ECL to the track momentum, the ECL shower 
shape, position matching between track and ECL clus- 
ter, the energy loss in the CDC, and the response of the 
ACC counters. Muons are identified based on their pen- 
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etration range and transverse scattering in the KLM de- 
tector. In the momentum region relevant to this analysis, 
charged leptons are identified with an efficiency of about 
90% while the probability to misidentify a pion as an 
electron (muon) is 0.25% (1.4%) [H,[2(|- Lepton tracks 
have to be associated with at least one SVD hit. In the 
laboratory frame, the momentum of the electron (muon) 
is required to be greater than 0.30 GeV/c (0.60 GeV/c). 
We also require the lepton momentum in the cm. frame 
to be less than 2.4 GeV/c to reject continuum. More 
stringent lepton requirements are imposed later in the 
analysis. 

For electron candidates we attempt bremsstrahlung re- 
covery by searching for photons within a cone of 3 degrees 
around the electron track. If such a photon is found, it 
is merged with the electron and the sum of the momenta 
is taken to be the lepton momentum. 



C. Background estimation 

Because we do not reconstruct the other B meson in 
the event, the B momentum is a priori unknown. How- 
ever, in the cm. frame, one can show that the B direction 
lies on a cone around the (£>*£)-axis [2l[ with an opening 
angle 2cos9b,D'£ defined by: 



cos ( 



'B,D*e = 



O TP * TP * 



z\Fb\\p1 



(15) 



In this expression, Eg is half of the cm. energy and \p£\ 
is sjE^f — m 2 B . The quantities E^, e , Pp, e and mD*i are 
calculated from the reconstructed D*£ system. 

This cosine is also a powerful discriminator between 
signal and background: Signal events should lie in the 
interval (—1,1), although - due to finite detector reso- 
lution - about 5% of the signal is reconstructed outside 
this interval. The background, on the other hand, does 
not have this restriction. 

The signal lies predominantly in the region defined by 
144 MeV/c 2 < Am < 147 MeV/c 2 and p e > 0.8 GeV/c 
(pi > 0.85 GeV/c) for electrons (muons). The region out- 
side these thresholds can be used to estimate the back- 
ground level. 

We therefore perform a fit to the three-dimensional 
(cos 9b,D'£, Am, pe) distributions. The cos9b,d*i range 
between —10 and 5 is divided into 30 bins. The Am (j>i) 
range is divided into five (two) bins, with bin boundaries 
at 141, 144, 147, 150, 153, 156 MeV (0.3, 0.8, 3.5 GeV/c for 
electrons and 0.6, 0.85, 3.5 GeV/c for muons). 

The background contained in the final sample has the 
following six components: 

1. D**: background from B -> D**£ + v decays with 
D** — > D*ir or D** — > Dir and from non-resonant 
B — > D*tt£ + v events, where the lepton has been 
correctly identified; 

2. correlated background: background from processes 
other than B — s- D**£ + i> decays in which the D* 



and the lepton originate from the same B meson, 
e.g, B° -> D*~t + v, t + ^ + vv\ 

3. uncorrelated background: the D* and the lepton 
come from different B mesons and the lepton is 
not from a B — >• D**£ + i> decay; 

4. fake lepton: the charged lepton candidate is a 
misidentified hadron while the D* candidate may 
or may not be correctly reconstructed; 

5. fake D*: the D* candidate is misreconstructed; the 
lepton candidate is identified correctly, but it is not 
from afi-» D**£ + v decay; 

6. continuum: background from e + e~ — » qq (q = 
u, d, s, c) processes. 

To model the D** component, which consists of a total 
of four resonant {D\, Dq, D[, DJ) and one non-resonant 
D*tt£i/ mode for both neutral and charged B decays, we 
reweight the branching ratios of each subcomponent to 
match the values reported by the Particle Data Group 0] ■ 
For the resonant parts, only products of branching ratios 
B(B D**£v) x B(D** -> DWtt) are available and 
consequently we reweight these products. The shape 
of the D** momentum distributions is also reweighted 
in 22 bins of q 2 to match the predictions of the LLSW 
model HHH. 

All of the background components are modeled by MC 
simulation except for continuum events; these are mod- 
eled by off-resonance data. For muon events, the shape 
of the fake lepton background is corrected by the ra- 
tio of the pion fake rate in the experimental data over 
the same quantity in the Monte Carlo, as measured us- 



ing K% -> 



decays. The lepton identification ef- 



ficiency is corrected by the ratio between experimental 
data and Monte Carlo in 27 — > ee//i/x events [H, l2(ij |. 
The (cos 6b, vi, Am, pi) distribution in the data is fit- 
ted using the TFractionFitter algorithm [24| within 
ROOT [251 ] . The fit is done separately in each of eight 
subsamples defined by the experiment range and the lep- 
ton type. The results are given in Table[J Figure |3] shows 
plots of the projections in co&6b,d*z for subsamplc B. 

In all fits, the continuum normalization is fixed to the 
on- to off-resonance luminosity ratio, corrected for the 
1/s dependence of the e + e~ — > qq cross section. In gen- 
eral, the normalizations obtained by the fit agree well 
with the MC expectations except for the D** component 
and the fake £ component, which are overestimated in 
the MC. After the background determination only can- 
didates satisfying the requirements —1 < cosOb.vi < 1, 
144 MeV/c 2 < Am < 147 MeV/c 2 and p e > 0.8 GeV/c 
(pi > 0.85 GeV/c) for electrons (muons) are considered 
for further analysis. 



D. Kinematic variables 

To calculate the four kinematic variables defined in 
Eq. Q]and Fig. [2]- W, cos 6^, cos6V and x ~ which char- 
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FIG. 3: Result of the fits to the (cos6 B D *o e , Am, pe) distributions in the e mode (left) and /i mode (right) of the sub-sample B. 
The bin boundaries are discussed in the text. The points with error bars are continuum-subtracted on-resonance data. Where 
not shown, the uncertainties are smaller than the black markers. The histograms are, top to bottom, the signal component, 
D** background, signal correlated background, uncorrelated background, fake i component and fake D* component. 




FIG. 4: Reconstruction of the B° direction. Refer to the text 
for details. 



acterize the B — > D*~l + v decay, we need to determine 
the B° rest frame. The B direction is already known 
to lie on a cone around the (-D*£)-axis with an opening 
angle 28b, we hi the cm. frame, Equation (|15l) . To ini- 



tially determine the B direction, we estimate the cm. 
frame momentum vector of the non-signal B meson by 
summing the momenta of the remaining particles in the 
event (p* ncl [HI) and choose the direction on the cone 
that minimizes the difference to —p* nci , Fig. |U 

To obtain p* nci , we exclude tracks passing far from the 
interaction point. The minimal requirements depend on 
the transverse momentum of the track, pt, and are set 
to dr > 20 cm (15 cm, 10 cm) or \dz\ > 100 cm (50 cm, 
20 cm) for a track p T < 250 MeV/c (p T < 500 MeV/c, 
Pt > 500 MeV/c). Track candidates that are compatible 
with a multiply reconstructed track generated by a low- 
momentum particle spiraling in the central drift chamber 
are also checked for and only one of the multiple tracks is 
considered. Unmatched clusters in the barrel region must 
have an energy greater than 50 MeV. For clusters in the 
forward (backward) region, the threshold is at 100 MeV 
(150 MeV). We then compute p[ nc \ (in the laboratory 
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A, Kit, e 


A, Kn, fj, 


B, Kn,e 


B, Kn,p 


Num. Candidates 


14802 


14203 


29217 


26894 


Signal events 


11609 ± 181 


11139 ± 190 


23029 ± 280 


21002 ± 258 


Signal (%) 


78.43 ± 1.22 


78.43 ± 1.34 


78.82 ± 0.96 


78.09 ± 0.96 


D" (%) 


5.63 ± 0.78 


4.02 ± 0.86 


4.32 ± 0.66 


3.90 ± 0.60 


Signal correlated (%) 


1.07 ± 0.17 


1.41 ± 0.25 


1.33 ± 0.16 


1.71 ± 0.19 


Uncorrelated (%) 


7.24 ± 0.35 


6.01 ± 0.40 


7.19 ± 0.31 


6.31 ± 0.29 


Fake £ (%) 


0.36 ± 0.17 


1.99 ± 0.34 


0.50 ± 0.17 


2.10 ± 0.23 


Fake D* (%) 


2.59 ± 0.12 


2.81 ± 0.13 


3.07 ± 0.11 


2.96 ± 0.10 


Continuum (%) 


4.68 ± 0.54 


5.32 ± 0.59 


4.77 ± 0.38 


4.93 ± 0.40 




C, Kn,e 


C, Kn, fj, 


D, Kn,e 


D, Kn, \x 


Num. Candidates 


22056 


20428 


15871 


14719 


Signal events 


17301 ± 240 


15513 ± 235 


12365 ± 189 


11469 ± 205 


Signal (%) 


78.44 ± 1.09 


75.94 ± 1.15 


77.91 ± 1.19 


77.92 ± 1.39 


D** (%) 


5.15 ± 0.71 


5.22 ± 0.71 


4.54 ± 0.72 


4.67 ± 0.86 


Signal correlated (%) 


1.56 ± 0.27 


2.07 ± 0.37 


2.01 ± 0.26 


2.73 ± 0.43 


t t i j i / fry \ 

Uncorrelated (%) 


6.35 ± 0.35 


6.01 ± 0.33 


7.33 ± 0.38 


6.30 ± 0.40 


Fake £ (%) 


0.75 ± 0.18 


2.26 ± 0.28 


0.30 ± 0.19 


1.68 ± 0.38 


Fake D* (%) 


2.86 ± 0.12 


2.69 ± 0.11 


2.89 ± 0.13 


2.80 ± 0.14 


Continuum (%) 


4.88 ± 0.45 


5.81 ± 0.51 


5.02 ± 0.53 


3.89 ± 0.49 



TABLE I: The signal yield and the signal and background fractions (given in %) for selected events passing the requirements 
| cos6 b ,d*i\ < 1, 144 MeV/c 2 < Am < 147 MeV/c 2 and pe > 0.8 GeV/c {p e > 0.85 GeV/c) for electron (muon) channels. 



frame) by summing the 3-momenta of the selected parti- 
cles, 

Pincl = ^2Pi > (16) 

i 

where the index i stands for all particles passing the con- 
ditions above, and transform this vector into the cm. 
frame. Note that we do not introduce any mass as- 
sumption for the charged particles. The energy com- 
ponent of pi nc i is determined by requiring -E* ncl to be 

^bcam = yfs/2. 

With the B° rest frame reconstructed in this way, the 
resolutions in the kinematic variables are found to be 
about 0.025, 0.049, 0.050 and 13.5° for w, cos 6^, cos0 v 
and x, respectively. 



IV. ANALYSIS BASED ON THE 
PARAMETERIZATION OF CAPRINI ET AL. 

A. Fit procedure 

Our main goal is to extract the following quantities: 
the product of the form factor normalization and \V c b\, 
F(l)\V c b\ (Eq. [TT]), and the three parameters p 2 , i?i(l) 
and i?2(l) that parameterize the form factors in the 
HQET framework (Eqs. ll"2"HL4")) . For this, we perform 



a binned \ 2 n t to the w, cos Qt, cos 9y and x distribu- 
tions over nearly the entire phase space. Instead of an 
unbinncd fit, we fit the one-dimensional projections of 
w, cos 6 1, cos 9y and \. This avoids the difficulty of pa- 
rameterizing the six background components and their 
correlations in four dimensions. In addition, the one- 
dimensional projections have sufficient statistics in each 
bin. However, this approach introduces bin-to-bin corre- 
lations that must be accounted for. 

The distributions in w, cos0£, cos#y and x are di- 
vided into ten bins of equal width. The kinematically 
allowed values of w are between 1 and 1.504, but we re- 
strict the fit range to values between 1 and 1.5. In each 
sub-sample, there are thus 40 bins to be used in the fit. 
In the following, we label these bins with a single index 
i, i = 1, . . . , 40. The bins i = 1, . . . , 10 correspond to 
the bins of the w distribution, i = 11, . . . , 20 to cos 0£, 
i = 21,..., 30 to cos6V, and i = 31,..., 40 to the x 
distribution. 

The number of produced events in the bin i, N% , is 
given by 

N$ h = N B oB(D*+ -S- D°ir + )B(D° ->■ if-f + )r fj0 r ! , 

(17) 

where N^o is the number of B° mesons in the data sam- 
ple, B(D*+ -> D°7r+) and B(D° -> R'-n+) arc the D* 
and D branching ratios into the final state under con- 
sideration t b o is the B° lifetime 0], and I\ is the 
width obtained by integrating Eq. [TU] in the kinematic 
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variable corresponding to i from the lower to the up- 
per bin boundary (the other kinematic variables are in- 
tegrated over their full range). This integration is numer- 
ical in the case of w and analytic for the other variables. 
The expected number of events N^ xp is related to iV* h as 
follows 



N 



40 

x p ^ ^ 



R tJ e 3 N 3 



bkgrd 



(18) 



Here, €j is the probability that an event generated in 
the bin j is reconstructed and passes all analysis cuts, 
and Rij is the detector response matrix, i.e., it gives 
the probability that an event generated in the bin j is 
observed in the bin i. Both quantities are calculated 
using MC simulation. J /V bkgrd is the number of expected 
background events, estimated as described in Sect. IIIlUI 
Next, we calculate the variance of of N^ xp . We con- 
sider the following contributions: the Poissonian uncer- 
tainty in iV* h ; fluctuations related to the efficiency, es- 
timated by a binomial distribution with N repetitions 
and known success probability e,; a similar contribution 
related to Rij using a multinomial distribution; and the 
uncertainty in the background contribution 7V l bkgrd . This 
yields the following expression for of , 



a, 



R%e*N? 



R 



2^(1 ^) (A rth )2 _ 



40 

E 

Rij(l — Rij) 2/i\rth\2 , r>2 e j0- 



iV dat 

eHN^Y+R 



Nmc 
a 2 (N* kstd ) . 



(19) 



The first term is the Poissonian uncertainty in Nf h . The 
second and third terms are the binomial and multino- 
mial uncertainties, respectively, related to the finite real 
data size, where iVdata (-^ da ta) ^ s the total number of 
decays (the number of reconstructed decays) into the fi- 
nal state under consideration (Kir, e or /i) in the real 
data. The quantities and Rij are calculated from a fi- 
nite signal MC sample (Nmc and iV^ c ); the correspond- 
ing uncertainties are estimated by the fourth and fifth 
terms. Finally, the last term is the background contri- 
bution cr 2 (7V l bksrd ), calculated as the sum of the different 
background component variances. For each background 
component defined in Section 1111 CI we estimate its con- 
tribution by linear error propagation of the results de- 
termined in Section IIII CI For continuum, we estimate 
the error in the on-resonance to off-resonance luminos- 
ity ratio to be 1.0%. These variances give the diagonal 
elements of the covariance matrix Cij . 

In each sub-sample we calculate the off-diagonal ele- 
ments of the covariance matrix Cy , i ^ j as Npij—NpiPj , 
where pij is the relative abundance of bin in the 

2-dimensional histogram obtained by plotting the kine- 
matic variables against each other, pi is the relative num- 
ber of entries in the 1-dimensional distribution, and N is 



the size of the sample. Covariances are calculated for 
the signal and the different background components in 
the MC samples, and added with appropriate normaliza- 
tions. 

The covariance matrix is inverted numerically within 
ROOT [25| and, labeling the electron and muon mode 
in each sub-sample with the index k, \ 2 functions are 
calculated, 



Xk 



(20) 



where N^ s is the number of events observed in bin i in 
data sample k. We sum these two functions in each sub- 
sample and minimize the global x 2 with MINUIT [261 ]. 

We have tested this fit procedure using generic MC 
data samples. All results are consistent with expectations 
and show no indication of bias. 



Investigation of the efficiency of low momentum 
tracks 



The tracking efficiency of the Belle experiment is repro- 
duced well by MC simulations for tracks with momenta 
above 200 MeV/c, which we refer to as "high momentum 
tracks". However, a significant portion of the momen- 
tum spectrum of the slow pions emitted in the D* decay 
lies below this boundary. For low momenta, the effects 
of interactions with the detector material such as mul- 
tiple scattering and energy loss become important and 
might lead to a deviation between data and MC in the 
reconstruction efficiency. 

We use one half of the reconstructed B —> D*lv sample 
to obtain corrections to the MC reconstruction efficiency 
in the low momentum range, measured using real data. 
The second half is used to perform the analysis with a 
statistically independent sample. The results of the back- 
ground estimation shown in Table U are those obtained 
in the samples used for the analysis. Both of the samples 
contain about 120,000 signal events. 

The sample used to investigate the efficiency of low 
momentum tracks is divided into a total of six bins in 
p Wa . The bin borders of the first five are 50 MeV/c, 100 
MeV/c, 125 MeV/c, 150 MeV/c, 175 MeV/c and 200 
MeV/c. The region beyond 200 MeV/c defines the sixth 
bin. By subtracting the background, we obtain an esti- 
mate of the signal in data and form the ratio with the 
signal in MC in each bin, = N^/N™ . 

The high momentum range is used as normalization, 
no efficiency correction is applied there. In the lower 
momentum bins we obtain the ratios p-n^.i = /;//max, 
which are identical to the ratio of reconstruction effi- 
ciencies in the bins i and the high momentum region, 
Ptr s ,i — &;Amax- We calculate this set of ratios for the 
electron and muon modes and form the weighted aver- 
age, separately for each of the four sub-samples. These 
values are applied as weights when filling the MC his- 
tograms to correct the reconstruction efficiency. 
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Most systematic uncertainties cancel out in the ratios 
p na i- Only the uncertainties in the various background 
components give a small systematic contribution to the 
uncertainty. 

This procedure assumes that the distribution of events 
in the p ng spectrum is identical for data and MC. How- 
ever, one of the aims of the analysis is to measure the 
form factor parameters that govern this distribution. 
Therefore, an iterative procedure is adopted: we calcu- 
late one set of corrections, apply them and perform the 
analysis to determine J"(l)|Vd,| and the form factor pa- 
rameters. We then calculate a new set of corrections us- 
ing these results and repeat the analysis. The changes of 
the parameters during this iterative procedure are small 
and vanish after the third iteration. We assign an addi- 
tional systematic uncertainty to our results based on the 
stability of the corrections against changes in the form 
factor parameters. As will be shown in Table IIII1 this is 
a negligibly small contribution. 

C. Results of the fits and investigation of the 
systematic uncertainties in the subsamples 

After applying all analysis cuts and subtracting back- 
grounds, a total of 123, 427 ± 636 signal events are used 
for the analysis, divided into a total of four experimental 
sub-samples as mentioned above. The result of the fit to 
these data is shown in Fig. [5] and Table |TTJ The x 2 per 
degree of freedom, x 2 / n -d.f, of all fits is good. Table [TT1 
also gives the x 2 probabilities or P- values, P x 2 . 

To estimate the systematic uncertainties in these 
results, we consider contributions from the following 
sources: uncertainties in the background component nor- 
malizations, uncertainty in the MC tracking efficiency, 
errors in the world average of B(D* + — > D°ir + ) and 
B(D° -> R-TT+) as well as in the B(B -> D**lv) compo- 
nents [7| , uncertainties in the shape of the w distribution 
of B -» D**iv events based on the LLSW model [Hj], 
uncertainties in the B° lifetime Q , and the uncertainties 
in the total number of B° mesons in the data sample. 

To calculate these systematic uncertainties, we con- 
sider 300 pseudo-experiments in which one of 15 param- 
eters is randomly varied, using a normal distribution. 
The entire analysis chain is repeated for every pseudo- 
experiment and new fit results are obtained, in total for 
4500 variations. One standard deviation in the pseudo- 
experiment fit results for a given parameter is used as the 
systematic uncertainty in this parameter. 

The parameters varied in the pseudo-experiments are 
as follows: 

1. The corrections on the tracking efficiencies for low 
momentum tracks are varied within their respec- 
tive uncertainties. To obtain the most conservative 
estimate, the uncertainties in different momentum 
bins are assumed to be fully correlated. Therefore, 
this component corresponds to a single parameter 
in the toy MC. 



2. The lepton identification efficiencies are varied 
within their respective uncertainties [l9l [20| . 

3. The normalization of the continuum background is 
not correlated with any of the other backgrounds, 
it is therefore varied individually within the uncer- 
tainty on the on- to off-resonance luminosity ratio, 
which is 1.0%. 

4. The normalizations of the remaining five back- 
ground components are varied within the uncer- 
tainties listed in Table HI while taking into account 
the correlations found in the background estima- 
tion described in Section lill CI 

5. Uncertainties in the composition of the D** com- 
ponent are accounted for by varying each of the 
components contributing to the D** background 
within the uncertainty reported by the Particle 
Data Group Q- For the resonant modes, this is 
the uncertainty in the branching fraction products 
B(B -> D**£v) x B(D** -> £>«7r); for the non- 
resonant mode, this is the uncertainty in B(B — > 
D*-kIv). 

6. In addition, the shape of the q 2 distributions of the 
D** components is varied according to the LLSW 
model [22J and the uncertainties on the model pa- 
rameters as determined in Ref. [23j |. 

7. The number of B°B° events is obtained from the 
product of the number of T(45) events in the 
sample with the branching fraction of T(45) to 

a B°B° pair. We vary the fraction / H / f 00 = 

B(T(AS) -> B+B-)/B(T(4S) -> B°B°) within its 
uncertainty [7(. This affects both the overall nor- 
malization and the background distributions. 

The uncertainties in the reconstruction of the high mo- 
mentum tracks, the branching ratios B(D* + — > D°ir + ) 
and B(D a —> K~tt + ), the number of T(45) events in 
the sample, and the B° lifetime affect only .F(l)|Vy, 
not the form factors. Therefore, their uncertainties are 
considered by analytical error propagation. 

D. Averaging the results of the subsamples 

To obtain the average of the four sub-samples, which 
have been measured independently, we use the algorithm 
applied by the Heavy Flavor Averaging Group [27| to ob- 
tain the world average for \V c b\ from semileptonic B de- 
cays. This algorithm combines both the statistical and 
the systematic uncertainties. The correlations of some of 
these errors between different samples is considered. For 
example, the uncertainty on the D* —> D°tt branching 
fraction will lead to a fully correlated systematic uncer- 
tainty in each B — > D*lv analysis. 

The average is obtained with the MINUIT package (2(| 
by using a \ 2 minimization. Here, N p gives the total 
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Sample A 


D° ->■ JTtt,! = e 


D° -> Kit, I = (j, 


total sample 


p 2 


1.248 ± 0.102 ± 0.022 


1.285 ± 0.114 ± 0.028 


1.259 ± 0.076 ± 0.019 


fli(l) 


1.317 ± 0.099 ± 0.041 


1.577 ± 0.131 ± 0.036 


1.436 ± 0.078 ± 0.030 


R 2 (l) 


0.804 ± 0.076 ± 0.017 


0.768 ± 0.093 ± 0.020 


0.795 ± 0.058 ± 0.015 


T(l)\V cb \ x 10 3 


34.8 ± 0.5 ± 1.2 


34.6 ± 0.6 ± 1.2 


34.7 ± 0.4 ± 1.2 


X 2 / n.d.f 


32.2 / 36.0 


31.6 / 36.0 


70.9 / 76.0 


P X 2 


0.651 


0.676 


0.643 


Sample B 


D° -»■ Kn,£ = e 


D° -+Ktt,£ = [j, 


total sample 


P 2 


1.169 ± 0.079 ± 0.011 


1.167 ± 0.088 ± 0.016 


1.168 ± 0.059 ± 0.011 


fli(l) 


1.411 ± 0.079 ± 0.026 


1.449 ± 0.090 ± 0.028 


1.427 ± 0.059 ± 0.022 


«2(1) 


0.902 ± 0.054 ± 0.011 


0.859 ± 0.061 ± 0.013 


0.882 ± 0.041 ± 0.010 


^(lJIVcil x 10 3 


34.4 ± 0.4 ± 1.1 


33.9 ± 0.4 ± 1.1 


34.2 ± 0.3 ± 1.1 


X 2 / n.d.f 


22.7 / 36.0 


36.5 / 36.0 


60.7 / 76.0 




0.958 


0.443 


0.900 


Sample C 


D° -^Kn,£ = e 


D° ^Kn,£ = fi 


total sample 


P 2 


1.226 ± 0.088 ± 0.011 


1.262 ± 0.101 ± 0.016 


1.239 ± 0.066 ± 0.011 


fli(l) 


1.363 ± 0.086 ± 0.026 


1.480 ± 0.107 ± 0.033 


1.411 ± 0.066 ± 0.023 


fla(l) 


0.891 ± 0.062 ± 0.012 


0.851 ± 0.076 ± 0.015 


0.876 ± 0.048 ± 0.012 


7-(l)|Vcfc| x 10 3 


34.4 ± 0.5 ± 1.1 


33.9 ± 0.5 ± 1.1 


34.2 ± 0.3 ± 1.1 


X 2 / n.d.f 


38.6 / 36.0 


38.2 / 36.0 


81.4 / 76.0 


P x2 


0.352 


0.370 


0.314 


Sample D 


D° -^Kn,£ = e 


D° -^Kn,£ = fi 


total sample 


P 2 


1.321 ± 0.102 ± 0.019 


1.174 ± 0.106 ± 0.020 


1.247 ± 0.073 ± 0.014 


fli(l) 


1.448 ± 0.109 ± 0.041 


1.230 ± 0.089 ± 0.031 


1.330 ± 0.069 ± 0.027 


J? 2 (l) 


0.791 ± 0.081 ± 0.019 


0.931 ± 0.071 ± 0.015 


0.864 ± 0.053 ± 0.014 


7-(l)|Vc6| X 10 3 


35.4 ± 0.6 ± 1.2 


35.7 ± 0.6 ± 1.2 


35.6 ± 0.4 ± 1.2 


X 2 / n.d.f 


25.1 / 36.0 


42.0 / 36.0 


70.1 / 76.0 




0.913 


0.226 


0.669 



TABLE II: The fit results for the four sub-samples. The first two columns show results obtained by investigating only the e or 
the fx channel, the third column is obtained by minimizing the sum of the \ 2 values calculated for each channel. Where given, 
the first error is statistical, and the second is systematic. 



number of fit parameters, in our case N p = 4. When 
calculating the average of N measurements of the four 
fit parameters pj = {F(l)\Vcb\, p 2 , Ri(l), i?2(l)}, a total 
of 4 x N values are available as inputs, which we la- 
bel as Vi — {Vi, . . . , V^n}. In general this number can 
be labeled as JVj. Each measurement Vi corresponds to 
one of the parameters pj , which defines a primitive map 
ir(i) : i — > p. The statistical covariance matrix of each 
measurement is known, as well as the correlation between 
the samples. The latter are zero in our case. This infor- 
mation allows one to construct a AN x 4iV-dimensional 
covariance matrix containing the statistical uncertainties 
and to obtain the statistical part of the x 2 to be mini- 
mized: 

xLt = E E ( V < ~ ^W) - V Aj) ) , (21) 

* 3 

where Vk indicates the average values of the fit parame- 
ters. 



The systematic uncertainties are implemented by as- 
suming Gaussian error distributions. The possible bias 
of input i with respect to the systematic source s can 
therefore be estimated as cr^s x r s , where r s is a normal 
distributed random number. Its square is by definition 
distributed according to a x 2 distribution with one degree 
of freedom, r 2 s ~ X 2 (l)- The parameters r s are floating 
fit parameters in MINUIT, and in each minimization step 
the sums 

Vi = V i + y £ j <j i , s r s ,Vi, (22) 



are evaluated. If a systematic uncertainty is associated 
with two different inputs i\ and they are both varied 
by the same fractional systematic uncertainty at the same 
time. The correlation between systematics is therefore 
included. 
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i-Resonance data 



Hated 



in 



°-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 °0 20 40 60 80 100 120 140 160 180 

cos e„ |x| [degrees] 



FIG. 5: Result of the fit of the four kinematic variables in the sub-sample B. The electron and muon modes are added in 
this plot. The points with error bars are continuum subtracted on-resonance data. Where not shown, the uncertainties are 
smaller than the black markers. The histograms are, top to bottom, the signal component, D** background, signal correlated 
background, uncorrelated background, fake £ component and fake D* component. 




FIG. 6: Plots of the result of the averaging procedure. Projections in J-"(1)|V C 6| vs. p 2 (top left), .F(l)|14&| vs. -Ri(l) (top 
middle), J-"(1)|T4j| vs. i?2(l) (top right), p 2 vs. -Ri(l) (bottom left), p 2 vs. R2(l) (bottom middle) and -Ri(l) vs. -R2(l) (bottom 
right) are shown. The red dot (solid line) shows the position (la ellipse) of the average, the blue rectangle (dashed line) the 
position (la ellipse) of the sub-sample A, the green triangle (dash-dotted line) the position (la ellipse) of the sub-sample B, 
the magenta diamond (dash-double dotted line) the position (la ellipse) of the sub-sample C and the cyan cross (dash-triple 
dotted line) the position (la ellipse) of the sub-sample D. 



The total \ 2 to be minimized takes the form, 



x 2 = 



3 3 



EE (ft- ^w) - fyi>) + E * 

S (23) 

and is minimized numerically. The number of degrees of 
freedom are calculated as 



any systematic uncertainties. The minimization is nu- 
merically stable and yields both the central values and 
the total uncertainties of the full four dimensional aver- 
age. 



n.d.f = (N t + N s ) - (N p + N s ) =Ni-N p , (24) 

floated parameters 

which is the same result one obtains in the case without 



Applying this procedure to the four results presented 
in Table [TT] yields the final result of this analysis. We 



12 





p 2 


fli(l) 


fla(l) 


T(l)\V cb \ x 10 3 


B(B° -> D*£v) [%} 


Value 


f.214 


1.401 


0.864 


34.6 


4.58 


Statistical Error 


0.034 


0.034 


0.024 


0.2 


0.03 


Systematic Error 


0.009 


0.018 


0.008 


1.0 


0.26 


Fast track efficiency 








-0.78 


-0.206 


Slow track efficiency 


+0.002 


+0.003 


-0.004 


-0.28 


-0.059 


Ptv s stability 


+0.00f 


-0.001 


+0.000 


-0.03 


-0.003 


LeptonlD 


+0.002 


+0.006 


-0.002 


-0.38 


-0.100 


Norm - D" 


+0.00f 


+0.001 


-0.001 


-0.03 


-0.008 


Norm - Signal Corr. 


+0.002 


-0.003 


+0.002 


+0.02 


+0.006 


Norm - Uncorr 


+0.002 


+0.008 


-0.003 


-0.02 


-0.001 


Norm - Fake I 


+0.003 


-0.003 


-0.001 


-0.01 


-0.003 


Norm - Fake D* 


+0.001 


-0.001 


+0.000 


+0.00 


+0.003 


Norm - Continuum 


+0.002 


+0.002 


-0.001 


+0.00 


-0.003 


D** composition 


+0.004 


+0.009 


-0.003 


-0.10 


-0.025 


D** shape 


+0.003 


+0.005 


-0.002 


-0.04 


-0.011 


iV(T(4S)) 








-0.24 


-0.063 


/+-//00 


+0.004 


-0.009 


+0.003 


+0.24 


+0.062 


B° life time 








-0.10 


-0.027 


B(D* D°7T S ) 








-0.13 


-0.034 


B(D° ->■ Kir) 








-0.22 


-0.059 



TABLE III: The breakup of the systematic uncertainty in the result of the fit to the full sample. The sign + (-) implies whether 
the fit result moves to larger (smaller) values, if the value of the corresponding systematic parameter is increased. 





Hi)\Vct>\ 


p 2 


fli(l) 


Jfe(l) 


■^(1)1^1 


1.000 


0.625 


-0.122 


-0.206 


p 2 




1.000 


0.575 


-0.872 


fll(l) 






1.000 


-0.697 


R 2 (l) 








1.000 



TABLE IV: The statistical correlation coefficients of the four 
parameters in the fit to the full sample. 



obtain 



•F(l)|Vc6| = 

P 2 = 

Mi) = 



(34.5 + 0.2 + 1.0) x 10" 
1.214 + 0.034 + 0.009, 
1.401 + 0.034 + 0.018, 
0.864 + 0.024 + 0.008, 



(25) 



with a x 2 /n.d.f = 14.3/12 (P x 2 = 0.282). This implies 
excellent agreement between the results, which can also 
be seen in the projections of the minimization, shown 
in Fig. [5] The corresponding branching fraction for the 
process B° — > D*~l + v is obtained from the integral of 
the differential decay width. We obtain 



B{B° -> D*~£ + v) = (4.56 ± 0.03 ± 0.26)%. 



(26) 



A breakdown of the systematic uncertainties is shown 
in Table [TTTT The statistical correlation coefficients of the 
result can be found in Table HVl 



MODEL-INDEPENDENT DETERMINATION 
OF HELICITY FUNCTIONS 



The angular distributions given in Eq. (|10[) are deter- 
mined by the kinematic properties of the decay. However, 
as discussed in section III C[ the expressions of the helic- 
ity amplitudes and thus the distribution in the variable 
w are based on the parameterization scheme proposed 
by Caprini, Lellouch and Neubert [3( . In this section, we 
extract the form factor shape of the longitudinal and the 
transverse components of Eq. [TO] through a fit to the w 
vs. cos 9y distribution. The binning is the same as in 
the fit approach described above. The contribution from 
events with w > 1.5 is fixed to the small values predicted 
by the results of the parameterized fit. 



A. Fit procedure 

From Eq. iflO)) we can obtain the double differential 
decay width dT/dw d cos 9y by integration over cos 0£ and 
X- 

If we define 



G 2 F {m B - m fl .) 2 mj), 

4 3 7T 3 



(27) 
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and 



B. Results 



,±± 



(w) = Vw 2 ~l(w + l) 2 h 2 Ai (w)\V cb \ 2 x 



1 — 2wr 



it 



w — 1 



(1 - r) 2 I ' V w + 1 

<7°V) = > A^(t U + l) 2 fc» i ( W )|V r c5 | 2 X 
iu — 1 



1 



1 - r 



(i-i? 2 H) 



(28) 



this equation becomes 

d 2 T{B+ -» £>*°l+^) 
c?w d cos 0v 



= F r [sin 2 0y (.g+ 4 
+2 cos 2 fly 5°°] 



(29) 



The quantities g correspond to the product of in- 
dependent parts of the different helicity combinations 
and kinematic factors. The one-dimensional distribution, 
as given in Eq. (fTT|l . depends only on the sum of these 
three combinations, 



dT{B+ -> D*°£+v t 
dw 



9 



9 



(30) 



The bin contents of the two-dimensional histogram in 
w vs. cos 6y can be obtained by integration of Eq. (|29|) 
over the corresponding bin area and considering the re- 
construction efficiencies and detector response as de- 
scribed in Eq. [T5] Each bin content can be given as 
the linear combination of two linearly independent parts. 
The integration of the angular distributions is performed 
analytically, the integration with respect to w defines a 
set of dimensionless parameters, 



G 



dw g 



kk 



(31) 



where Wj — {w\, W2, ■ ■ ■ , nn} = {1, 1.05, . . . , 1.5} are the 
bin boundaries of the 10 bins in w. In addition, wc define 



->++ 



G 



C++ + G: 



9 



9 



oo 



9 = 9 
r<L _ /^<oo 
u i — u i . 

For the w vs. cos#y distribution we the x 2 function 



and 



xt 



EE 

»=1 3=1 



a r2D,obs a r2Z),cxp 
k,ij k,ij 



Nf 



(32) 



for each channel k, which depends only on the parame- 
ters Gf and G\ . Here N° hs gives the number of events 
observed in on-resonance data, N e * p the number of ex- 
pected events, as defined in Eq. 1181 and <7atcx P the un- 
certainty in the expected number of events, as given in 
Eq. (fl9| . Again we form the sum of % 2 functions for each 
channels and minimize this expression numerically using 
MINUIT [26j . 

We have tested this fit procedure using generic MC 
data samples. All results are consistent with expectations 
and show no indication of bias. 



We investigated the largest of the four subsamples - 
corresponding to about 274 fb _1 of data recorded on the 
T(45) resonance - to extract the helicity shapes. Tables 
IVl and I VII give the results of the fits, where the system- 
atic errors quoted in these tables originate from the same 
sources given in the breakdown in Table IIII1 Many of the 
dominant systematic uncertainties, such as the track re- 
construction or the lepton ID uncertainty, are correlated 
between different bins. The \ 2 °f the fit is statistically 
consistent with the number of degrees of freedom, we ob- 
tain x 2 /n.d.f = 175.8/179, P x 2 = 55.4%. The results are 
shown in Fig. [7] and Tables Ivll VII compared to the values 
obtained using the parameters given in Section HV Dl 

Finally we compute a numerical measure of the agree- 
ment between the result obtained from this two dimen- 
sional fit with the functions predicted by the parameter- 
ization of Caprini, Lellouch and Neubert Q • Comparing 
with the values given in Table Hi] is naturally problem- 
atic, since the underlying sample is identical in both fits. 
Therefore we compare the extracted shapes to the pa- 
rameterization using the world average reported by the 
Heavy Flavor Averaging Group 27] in the summer of 
2010 and form a simple x 2 , 



X 



LLP 



par am 



Gf) (G 



i ) 

stat / 



(G 



p. ram ^ .fit j 



(33) 

where (jP aiam arc the values obtained using the param- 
eterization by Caprini, Lellouch and Neubert, Gf* are 
the values obtained by the unparameterized fit described 
above and G s t a t is the covariance matrix between the 
fit parameters, which is also a result of the numeri- 
cal minimization. Evaluating this expression yields a 
X 2 /n.d.f = 29.3/20, P x 2 = 8.3%. This implies satis- 
factory agreement between the extracted shapes of the 
longitudinal and transverse components and the world 
average parameters. 



VI. SUMMARY AND DISCUSSION 

Using 711 fb _1 of data collected by the Belle exper- 
iment, we have analyzed approximately 120,000 B° — > 
D*~l + vi decays. A fit to four kinematic variables 
fully characterizing these decays yields measurements of 
the product of the form factor normalization and IV^I, 
■7-"(l)|Vd>|, ancl °f the parameters p 2 , i?i(l) and i?2(l) 
that enter the HQET form factor parameterization of 
this decay. We obtain: 

T(l)\V cb \ = (34.6 ± 0.2 ± 1.0) x 10" 3 , 

p 2 = 1.214 ±0.034 ±0.009, 

Rx{l) = 1.401 ±0.034 ±0.018, 

R 2 {1) = 0.864 ± 0.024 ±0.008, 

B(B° D*-£ + u e ) = (4.58 ±0.03 ±0.26)%. (34) 
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TABLE V: Results obtained for Gj (dimensionless) , compared to the central values obtained from the parameterized fit. 
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TABLE VI: Results obtained for Gf (dimensionless), compared to the central values obtained from the parameterized fit. 
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FIG. 7: Results of the fit of the helicity functions (red crosses) compared to the prediction obtained by using the parameters 
obtained by using the parameterization prescription by Caprini, Lellouch and Neubert (solid black line). The left plot shows 
the results for Gf , the right one for Gj . Only the statistical errors are shown. 



For all these measurements, the first error is the statis- 
tical uncertainty and the second is the systematic un- 
certainty. Using a recent lattice QCD result, J-{1) — 
0.921 ± 0.013 ± 0.020 [loj], we obtain the following value 

Of |Vc6|, 

I V cb \ =37.5 ±0.2 ±1.1 ±1.0, (35) 

where the third error is due to the theoretical uncer- 
tainty on J-"(l). Our results (j3~4"]) are compatible with the 
recent measurements of these quantities by the BaBar 
experiment 12814301 as well as with results reported by 
the ALEPH [31(, CLEO H, DELPHI [Ht H and 
OPAL [35| experiments. This paper supersedes our pre- 
vious result [36|, based on a subset of the data used in 
this analysis. The results presented here give the most 
precise determination of the form factor parameters and 
J"(l)|F c6 | to date. 

A direct, mo del- independent determination of the form 
factor shapes has also been carried out and is in good 
agreement with the HQET-based form factor parameter- 
ization by Caprini et al. Q. 
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